Rydberg-London Potential for Diatomic Molecules and Unbonded Atom Pairs 
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We propose and test a pair potential that is accurate at all relevant distances and simple enough for 
use in large-scale computer simulations. A combination of the Rydberg potential from spectroscopy 
and the London inverse-sixth-power energy, the proposed form fits spectroscopically determined 
potentials better than the Morse, Varnshi, and Hulburt-Hirschfelder potentials and much better 
than the Lennard- Jones and harmonic potentials. At long distances, it goes smoothly to the correct 
London force appropriate for gases and preserves van der Waals's "continuity of the gas and liquid 
states," which is routinely violated by coefficients assigned to the Lennard- Jones 6-12 form. 



There are at least three classes of interatomic poten- 
tials. Commercial codes use the Lennard- Jones 
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forms, which are accurate near the minimum at r = ro- 
But this first class of potentials may be too simple for 
complex materials away from from equilibrium. 
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(J3 = (nr 2 - l)/(2r 2 )), and Hulburt-Hirschfelder 

V HH (r) = V M (r) + \V(r )\ c'k 3 x 3 e" 2ra (l + Kb'x) (5) 

(for b', c', see ref.0) potentials represent a second class 
of potentials Q accurate over a wider range of distances. 

Quantum chemists 0, El 01 El El have derived a third 
class that reproduce spectroscopic and thermodynamic 
data with impressive fidelity. But the potentials of this 
class involve many parameters and may be too cumber- 
some for use in large-scale simulations. 

We propose and test a form 



V(r) = ae~ br (1 - cr) 
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that is nearly as accurate as the class-3 potentials but 
simpler than many class-2 potentials. It is a combina- 
tion of the Rydberg formula used in spectroscopy and 
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FIG. 1: The hybrid form V (with a = 53.8 eV, b = 2.99 A" 1 , 
c = 2.453 A" 1 , d = C 6 = 3.884 eVA 6 , and e = 47.6 A 12 ) 
(solid, red) fits the RKR spectral points for the ground state 
of molecular hydrogen (pluses, blue) and gives the correct 
London tail for r > 3 A. The Lennard- Jones Vlj (dashes, 
green) and harmonic Vh (dots, magenta) forms fit only near 
the minimum. 



the London formula for pairs of atoms. In Eq.@), the 
terms involving a. b, and c were proposed by Rydberg 01 
to incorporate spectroscopic data, but were largely ig- 
nored until recently 0, The constant d = C§ is 
the coefficient of the London tail. The new term er~ 6 
cures the London singularity. As r — > 0, V(r) — > a, 
finite; as r — ► oo, V(r) approaches the London term, 
V(r) — > —d/r 6 = —Cs/r 6 . In a perturbative analy- 
sis |l3j . the a,b,c terms arise in first order, and the d 
term in second order. 

To test whether the hybrid V(r) can represent covalent 
bonds far from equilibrium, we used Gnuplot 0] to fit 
Eq.(JSJ) to empirical potentials for molecular H2, N2, and 
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FIG. 2: The hybrid form 1/ (Eq.JgJ as in Fig.[T] solid, red) fits 
the RKR spectral points for the ground state of molecular hy- 
drogen (pluses, blue). The Morse Vm (Eq.©, dashes, green), 
Varnshi Vv (Eq.(|IJl, dots, magenta), and Hulburt-Hirschfelder 
Vhh (Eq.@, dot-dash, cyan) are too low for 1.5 < r < 3 A. 
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FIG. 4: For the ground state of molecular nitrogen, the hybrid 
form V (Eq.@, as in Fig.|3] solid, red) fits the RKR spectral 
points (pluses, blue) and the FO points (x's, black). Vm, Vv, 
& Vhh as in Fig.H 
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FIG. 3: For the ground state of molecular nitrogen, the pro- 
posed potential V of Eq.© (with a = 3.25 keV, b = 4.30 A" 1 , 
c = 1.1906 A -1 , d = 14.1 eVA 6 , and e = 27.1 A 12 ) (solid, 
red) fits the RKR spectral points (pluses, blue) and the FO 
points (x's, cyan). Vlj and Vh as in Fig. Q 



O2 obtained from spectroscopic data |3 0, 0] by the 
RKR (Rydberg [ljj, Klein 0, Rees 0) method, set- 
ting d equal to the London values. Figure shows that 
the hybrid potential of Eq.© (solid, red) goes through 
the RKR points for H2 (pluses, blue) from 0.5 to 4 A. Fit- 
ted to the minimum, the harmonic potential (dashes, 
green) and the Lennard-Jones potential (dashes, blue) 



> 

CD 



1 






1 ; 1 1 1/ 1 

\v LJ jv H 




-1 












\ ', ;' /' / 








+ 1 y AS 












-2 












\ \ 












-3 












'■ \ 1 / ■ jf 








'■\ \ /'/ 




-4 




\ \ '. / / 








\ \ o 2 












-5 










0.8 



1.2 



1.4 

r(A) 



1.8 



FIG. 5: For the ground state of molecular oxygen, the poten- 
tial V of Eq.© (with a = 3.61 keV, b = 4.48 A -1 , c = 1.05 
A -1 , d = 16.08 eV A 6 , and e = 58.4 A 12 ) (solid, red) fits the 
RKR points (pluses, blue). Vlj and Vh as in Figs.0&|3] 



are accurate only near the minimum at 0.74 A. Fig- 
ures |21 & El show that the hy brid potential V{r) fits the 
N 2 and 2 RKR points [lH E3 and the N 2 first-order 
(FO) points between 1 and 1.8 A. The Lennard- 
Jones and harmonic potentials of Eqs.imS) n * only near 
the minima. 

How does the hybrid form compare with the class-2 
potentials of Eqs.GH5|l? Figures [3 H & © show that 
for 1.4 < r < 2 A the hybrid V(r) is closer than 
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FIG. 6: For the ground state of molecular oxygen, the hybrid 
form V (Eq.©, as in Fig.0 solid, red) fits the RKR spectral 
points (pluses, blue). Vm, Vv, & Vhh as in Fig. 
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FIG. 8: Positive potential V (Eq.JjJ as in Fig. 0, Ar 2 Aziz 
potential (pluses, blue), L-J form Vlj (Eq.Q as in Fig. 0. 
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FIG. 7: The potential V (Eq.@ for the Ar-Ar ground state 
with a = 1720 eV, b = 2.6920 A -1 , c = 0.2631 A" 1 , d = 
37.943 eVA 6 , and e = 177588 A 12 ) (solid, red) fits the Ar 2 
Aziz potential (pluses, blue) with the correct London tail. 
When matched at the minimum, the Lennard- Jones form Vlj 
(Eq.JTJ with r = 3.757 A and V(r ) = -0.01234 eV) (dashes, 
green) is too low for r > 4 A. 
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FIG. 9: The potential V (Eq.® for the Kr-Kr ground state 
with a = 2499 eV, b = 2.5249 A" 1 , c = 0.2466 A~\ 
d = 78.214 eVA 6 , and e = 199064 A 12 ) (solid, red) fits 
the Kr 2 Aziz points (pluses, blue) with the correct London 
tail. When matched at the minimum, the Lennard- Jones form 
Vlj (Eq.lJTJ with r = 4.008 A and V(r ) = -0.017338 eV) 
(dashes, green) is too low for r > 4 A. 



|SJ) to the RKR points. A useful estimate of how well a 
particular potential Vp(r) fits N data points Vd(Ti) is 
the dimcnsionlcss error 
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For H2, N2, and O2, the average error 5 was 59.7 for Vh, 



49.3 for Vlj, 0.037 for V M , 0.031 for V v , 0.021 for Vhh, 
and 0.0044 for the hybrid V. The hybrid form V is five 
times more accurate than the class-2 potentials (PHSJ and 
four orders of magnitude more accurate than the class- 1 
potentials (Q0). 

Can V(r) also represent weak noncovalent bonds? Us- 
ing Gnuplot, we fitted a, 6, c, and e in Eq.JSJ to Aziz's 
accurate HFDID1 potential for Ar 2 @ and HFD-B po- 
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FIG. 10: Positive potential V (Eq.@ as in Fig.®, Kr 2 Aziz 
points (pluses, blue), L-J form Vlj (Eq.Q as in Fig.®. 

tential for Kr 2 [13] and set d equal to their London-tail 
coefficients. Figures 1711101 show that for 2 < r < 7 A, the 
hybrid form © (solid, red) fits the Aziz potentials for 
both Ar2 and for Ki~2 (pluses, blue). The Lennard- Jones 
Vlj curves £Q| (dashes, green) matched at the minima 
are too deep for r > 4.5 A (Figs. H&IHJ) and too hard 
for r < 3 A (Figs. 131 fcEnj). The potential V(r) of Eq.© 
represents weak noncovalent bonds better than Vlj (and 
V H ). 

Does it matter that Vlj{t) fails to fit the Ar-Ar and 
Kr-Kr interactions? To find out, we used V and Vlj to 
compute the dimcnsionless second virial coefficient B 2 /tq 
of Ar and Kr at room temperature (kT — 0.025 eV). Here 
ro is the minimum of the potential, and the second virial 
coefficient B 2 is the integral over all space 

S 2 (T) = -I Jd\ (e-W-l) (8) 

in which f3 = l/(kT). The hybrid potential V fitted to 
the curves of Figs. 171- HUI gives B 2 /r% = -0.499 for Ar 
and —1.35 for Kr, which respectively differ from the ex- 
perimental values of —0.552 and —1.41 by 9.6% and 
4.1%. The Lennard- Jones potential Vlj fitted to ro and 
V(r ) gives B 2 /r% = -0.899 for Ar and -1.92 for Kr (er- 
rors of 63% and 37%); it also wags a long-range tail with 
London coefficients 2 |V"(ro)| r$ that are too large by 83% 
for Ar and by 84% for Kr. If the parameters ro and V(r<j) 
in Eq.Q for Vlj are chosen to give the correct second 
virial coefficient B 2 jr\ for a range of temperatures [23j |. 
then V(r ) is too shallow by 16 % for Ar and 15% for 
Kr, and the London coefficients of the long-range tail are 
too large by 70% for Ar and by 64% for Kr. With only 
two parameters, Lennard-Jones fits are procrustcan. 
What about additivity (HHHH^)? When three or 



more atoms interact, their potential energy is not the sum 
of the three (or more) pair potentials. Is the accuracy 
of the hybrid form important in the liquid phase where 
additivity is only approximate? 

To test whether the lack of complete additivity in the 
liquid phase obscures the advantages of the hybrid form 
V over the Lennard-Jones potential Vlj, we used both to 
compute the heats of vaporization A vap i? of Ar and Kr 
at their boiling points at atmospheric pressure. In our 
Monte Carlo simulations, we imposed periodic boundary 
conditions to reduce finite-size effects. Our Monte Carlo 
code is available at bio.phy s.unm.edu/latentHeat| We 
used it to compute the potential energy U per atom in 
the liquid and gas phases. The latent heat of vaporization 
A vap H is the difference between the potential energies 
Ug&s arid C^iq U j c j plus the work done in expanding by 
AV against the pressure p of the atmosphere, A vap iJ = 
t/gas - C/i iqmd + pAV. 

The hybrid form V fitted to the curves of Figs. 171- 1TUI 
gave A vapj ff = 0.0694 eV (per atom) for Ar and 0.0982 
cV for Kr, which differ from the experimental values 
of 0.0666 and 0.0941 eV by 4.2% and 4.4%. In equiva- 
lent Monte Carlo simulations, the Lennard-Jones poten- 
tial V L j fitted to r and V(r ) gave and A vap H = 0.0787 
cV for Ar and 0.111 eV for Kr (errors of 18% and 18%). 
So the errors due to a lack of additivity are of the order of 
4%, while those due to the Lennard-Jones potential are 
about 18%. Even in the liquid phase, limited additivity 
is less of a problem than the defects of the Lennard-Jones 
potential. 

For a wide class of atom pairs, the hybrid form © 
can reliably represent the best spectroscopically deter- 
mined potentials over all relevant distance scales. It also 
yields accurate second virial coefficients and heats of va- 
porization. Its simplicity recommends it as a teaching 
tool and as a practical form for computation. Given the 
differences between it and the Lennard-Jones, harmonic, 
Morse, Varnshi, and Hulburt-Hirschfcldcr potentials, it 
would be worthwhile to examine the consequences of 
these differences in Monte Carlo searches for low-energy 
states of biomolecules and in numerical simulations of 
phase transitions and reactions far from equilibrium. 
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